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N. DEMNI I 

Abstract. We supply two different descriptions of the pushing process 
driving the reflected Brownian motion in Weyl chambers, when the latter 
domains are simplexes. The flrst one shows that a simple root lies in one 
and only one orbit if and only if the pushing process in the direction of 
that simple root increases as the sum of all the Brownian local times in 
the directions of the orbit's positive elements. The last one shows that 
the pushing process may be written as the sum of an inward normal 
vector at the chamber's boundary and an inward normal vector at the 
origin, yielding a kind of a multivoque stochastic differential equation for 
the reflected process. We finally give a particles system interpretation of 
the reflected process and construct a multidimensional skew Brownian 
motion. 

1. OVERVIEW 

By 'the reflected Brownian motion' (shorthand RBM), it is often meant 
the absolute value of a real Brownian motion. This process has gained much 
fame since it is not an Ito's semimartingale and due to the celebrated Levy's 
representation in relation with Skorohod's problem ([lOj). More precisely, 
if 5 is a standard real Brownian motion, then there exist a standard real 
Brownian motion /3 and an increasing process L°(S) such that 

(1) \Bt\=l3t + L^t{B). 

The process LP is known as the local time at of since the following 
representation holds (Ch.VI in |10| ) 

and it gives the explicit Doob-Meyer's decomposition of the submartingale 
\B\. It also shows that L^{B) is adapted with respect to the natural filtration 
^1^1 of \B\. Besides, the support of dLl^{B) is contained in the set {t^Bt = 
0} and the identity in law of trajectories due to Paul Levy holds (^lOj) 

{\B\,L\B)) = {S-B,S), St := sup Bs. 

0<s<t 
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Motivated by applications to queueing theory and constructions of sta- 
ble subordinators, various probabilists were interested in multidimensional 
analogs of the RBM. Such processes behave as multidimensional Brownian 
motions in the interior of their state spaces and perform reflections at the 
boundaries. Examples of state spaces already appeared in literature include 
wedges (|13]). two dimensional cones 

{z = re'^ G C, |6l| < a, < a < tt/2} 

and more generally simplexes 

n 

{x G M", X = XiVi, Xi > 0} 
1=1 

with respect to some basis {vi)i so that a side, say Aj, 1 < z < n is specified 
by Xj = (HI), domains delimited by a set of hyperplanes (pjj), convex 
polyhedral domains: 

{x G M", {ni,x) > bi, I < i < d} 

where (nj)i<j<rf are n-dimensional vectors and (6j)i<j<rf is a d-dimensional 
vector ([!]) and in particular orthants ([I2j)- Most of the papers cited above 
were mainly intended to decide whether the RBM is a semimartingale or 
not and whether the uniqueness in law holds in the affirmative case. 
Recently, O. Chybiryakov defined and studied a RBM valued in a class 
of cones that arises in Lie algebra theory and more precisely associated 
with root systems: Weyl chambers ([2], [7]). From their very definition (see 
below), these domains are convex simple polyhedrons for which (nj)j are 
the so-called simple roots and bi = for all i, are simplexes for most of 
the root systems, in particular are orthants for the so-called orthogonal root 
system and are wedges of angles vr/m, m € {2,3 . . . } for dihedral systems 
(12J)- When it is not a simplex, the Weyl chamber may be a half plane 
(root systems of type A2), a setting that was of considerable interest to 
various mathematicians (see p. 406 in [13] for references). Assuming that 
Weyl chambers are simplexes, the author actually derived a Tanaka's type 
SDE and a representation of its drift, the pushing process, as an additive 
functional. With regard to the wealthy and regular geometrical structure 
of Weyl chambers and motivated by applications to queueing theory, we 
thought it is natural to answer the following questions: 

• Given a simple root, is it possible to write the coordinate of the 
pushing process in this direction as a sum of local times of some real 
semimartingales? 

• Is there any relation to multivoque SDKs, a tool that showed to be 
powerful in stochastic analysis with reflections ([I])? 

The answer to the first question is positive if and only if no other simple root 
lies in the orbit of the given simple root, while there is no perfect relation 
to multivoque stochastic analysis and this fact is not surprising since the 
reflections are not normal to the boundary. Nevertheless, we show that the 
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additive functional decomposes as the sum of an inward normal vector at 
the boundary of the Weyl chamber and of an inward normal vector at the 
vertex of the simplex. As a by product, the RBM satisfies a multivoque- 
type SDE with a repelling to infinity which agrees with the transience of the 
reflected process. The paper is closed by further developments: we express 
the pushing process by means of local times of the RBM and we give a 
construction of a multidimensional the skew Brwonian motion (^). For 
sake of completeness, we collect below some facts on root systems and fix 
some notations we use later on (Ch.VI in [5j, Ch.I, II in [7j). Once we do, 
we state Chybiryakov's non published result and we review its proof. This 
is two-fold: we correct some minor errata and we discuss the case of the root 
system of type A for which C is not a simplex of V. 

2. RBM IN Weyl chambers: Reminder 

2.1. Facts on root systems. Let {V,{,)) be a n-dimensional Euclidean 
space and define the reflection 

(ct v) 

a^{v) = v -2j-^a, a€V\{0} 

{a, a) 

with respect to the hyperplane Ha = O"*" := {x, {a,x) = 0}. Then, a root 
system i? is a collection {a £ ^\{0}} such that (Ta{R) = R for all a £ R. It 
is said to be reduced if Man i? = {ita}. A trivial example is provided by the 
orthogonal root system R = {±€{,1 < i < n}, where {ei)i is the canonical 
basis of V. It can be shown that there exists a basis S of span(i?) C V, called 
simple system, such that a root a £ Ris either a positive or a negative linear 
combination of elements of S. Thus R splits into a positive and a negative 
systems R+,R- respectively. With the help of elements of S, known as 
simple roots, one defines the so-called dosed Weyl chamber 

C:={x£ V, {a,x) > Ofor aha G 5}. 

This is a cone, not necessarily a simplex, and specializes to the positive 
orthant of M" when R is the orthogonal root system. Moreover, C enjoys the 
following important property: let W be the group spanned by {aa, a £ R}: 
the reflections group (it is a subgroup of the orthogonal group 0{n)). Then, 
C is a fundamental domain with respect to the action oiWonV{x^ wx), 
in the sense that for any x £ V, there exist w £ W and a unique y £ C such 
that x = wy. Thus, the space V decomposes as 

and the projection t: : V ^ C is well defined (note that t: = w* 

if X G wC for some w £ W), where is the orbits space. When S = 
{si, . . . , Sn} is a basis of V then C is a simplex of vertex Oy since 

n 

C = {x £V, X = '^{si,x)^i, {si,x) > 0, 1 < i < n}, 

i=l 
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where is the dual basis of S defined by: 

In the last decade, an active probabilistic research related to root systems 
has emerged after the introduction of the so-called Dunkl processes and 
their W^-invariant parts ([3]). The latter processes are diffusions valued in 
C and specialize to the RBM in Weyl chambers defined as vr(0), where Q 
is a y- valued Brownian motion. 

2.2. On Chybiryakov's result and review of its proof. 

Theorem 1. Assume that S = {si, . . . , Sn} is a basis ofV. Then there exist 
a V -valued Brownian motion B and a C -valued continuous process denoted 
LO(e) such that 

(2) 7r(et) = ^(eo) + i?t + LO(e), 

and for any 1 < i < n, {si, L^{@)) is an increasing process satisfying 

(3) ri|(,^,,(e,))^o}(s*,«!i°(0)) =0. 
Jo 

More precisely, if A is the matrix whose rows are the simple roots, then 

AL\e) = Y 

where Y = {Y)i<i<n with 

/g. g.\ r 

(4) Y = lim — LLJ_ \ / 1, , tn, \^ r- r^\ds, 1 < i < n. 

Remarks. 1/ Following the terminology used in literature, we refer to L^{@) 
as the pushing process of tt{@). Besides, we shall use the notation for 
both the pushing process and local times of real semimartingales, neverthe- 
less the argument between brackets indicates which process we refer to. 
2/ We point the readers interested in Chybiryakov's dissertation to the fact 
that the author assumed that all roots have equal length \f2. That is why a 
rescaling factor showed up in each coordinate Yi . 

3/ From results derived in [4J, one easily sees that the RBM in Weyl cham- 
bers is a SRBM in the sense of Williams and Dai, and uniqueness in law 
holds for the SDE 

4/ When S is a basis of V , the root system R is said to have full rank. 
However, this assumption is not always valid as it is not for the root system 
of type A defined by ([7\ p. 41) 

R = {±(ej — Cj), 1 < i < j < n}, S = {ci — e^+i.l < i < n — 1}. 

As the reader may easily check, the simple system S is no more a basis of 
V = M", rather spans the hyperplane orthogonal to the vector (1, . . . , 1) so 
that Aisanxn — 1 matrix and does not have a left inverse. Nevertheless, 
the same lines of Chybiryakov's proof lead to 

An{@t) = A7r(eo) + ABt + Y, 
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where Y = (yi)i<i<n_i is defined in the same way as in 

Review of the Proof. O. Chybiryakov used a smooth approximation of the 
projection vr together with Ito's formula to derive ([2]). But there is an errata 
in the expression of vr displayed in [3] p. 170. In fact, given x £ V, there are 
as many elements w as the cardinality of the isotropy subgroup of x such 
that w*x G C. Thus, w is unique if x G wC (since the isotropy group is 
trivial, [7J p. 22) and ir = w* and more generally one writes 

''^''^ = ^{weW,wx = x} ^ ^'^{-e-c^} 

which is no more differentiable (since the cardinality of the isotropy group 
is an integer- valued function). However, when one deals with a multidi- 
mensional Brownian motion which may hit one and only one hyperplane 
([6], p. 255), the isotropy group contains two elements: the identity and the 
reflection orthogonal to that hyperplane thereby 

at any time t. In this way, Chybiryakov's proof remains valid since {t, w*Qt £ 
dC} has zero Lebesgue measure by the occupation density formula ([2] p. 68), 
and it follows exactly the lines of the one written in [3] p. 172-175 with zero 
multiplicities. 

3. The pushing process as a sum of local times 

Let Si G S, then it is natural to ask whether Yi may be written as a 
sum of the local times L^{{'wsi,Q)),'w G W. This claim was shown to be 
true when R = B2 (p] p. 70) but this case is not illustrative since B2 is the 
most elementary example of a non orthogonal root system. Here we supply 
the appropriate necessary and sufficient condition ensuring the validity of 
this claim, thereby recovering Chybiryakov's result and covering other root 
systems like even dihedral ones ([5] p. 144). 

Proposition 1. Let e > 0, then Si is the only simple root in its orbit i?* = 
{wSijW G W} if and only if 

^wew{Q < {wsi,x) < e,x G wC} = U„gj^»nfi+{Ka> < e} 

for any x £ V. 

Proof. Assume there exists a simple root Sj,j 7^ i lying in R^. Then, 
Sj~ 7^ defines one of the hyperplanes delimiting C. Since every chamber 
wC is visited only once as w runs over W (counting intersections), then it 
is impossible to have < {sj,x) < e,x £ C if one already has < < 
e, X G C since for, C will be visited twice. Hence the former set is strictly 
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included into the latter one. Now, assume Si is the only root in R^, then the 
boundary of each chamber wC contains exactly one hyperplane Ha = a"*" 
for a = wsi. Since one visits any chamber wC exactly once as w ranges W, 
then one comes close to all hyperplanes for a = wSi G i?*, at a distance 
at most equal to e. This proves that 

U^gVK{0 < {wsi,x) <e,xe wC} = Ua(=^!{0 < {a,x) < e}. 

Since —a G i?*, the Proposition is proved. ■ 

Corollary 1. If Si is the only simple root in its orbit, then 

Y, = {s„L'{Q))= L'{{a,e)). 

Proof, for any a £ K^CiR^ and any e > 0, define Aa{x) := {|(a,x)| < e}. 
Then, one easily shows using the occupation density's formula (Ch.VI \1Q\) 
that 

/ . / ^{o<(s„w*es)<t,esewC}^^ = / ^u^Gw{o<{s^,w*e,)<6,esewC}ds 

..^w^O ' Jo 







{A^{es)\{a,0,)=O)}ds 



^^o.emnR+{Ac,{es)}ds. 
Now, one splits the lastly written indicator function into 

and all the possible intersections of sets Aa{Qs), s being fixed. Since (a, a) = 
{si,Si) for any a £ and since {a,Q) is a local martingale with brackets 
process {{a,a)t)t, then 

lim^^l^„(e.)d. = LO((«,G)) 

by Corollary 1.9 Ch.VI in [lOj. Now, each intersection of sets Aa{Q),a £ 
R^ n i?+ contains two (at least) non necessarily independent Brownian mo- 
tions (up to scaling factors). Thus, one constructs a planar Brownian motion 
whose length is less than ce for some positive constant c. As a matter fact, 
the indicator function of any intersection of sets Aa{Qs),0! G K^Ci R+ is less 
then 

l{0<Zs<e} 

where Z is a Bessel process of dimension 2. Since the local time of the latter 
process is almost surely zero, then all integrals of the indicator function of 
any intersection of sets Aa{Q),a G i?* n R^ vanish after rescaling by 1/e 
and letting e ^ 0. ■ 
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Remark. When two simple roots are conjugated, one can only claim 

U«;Giy{0 < (wsijx) < e,x £ wC} = U«,gvy{— e ^ {wsi,x) < e,x £ w(CUSjC)}. 

An illustrative example is provided by the odd dihedral root system /2(3) f|5] 
p. 144)' there are six two-dimensional roots represented as 

{±e-*''/2e^'"''/^ m= 1,2,3} 

and there is only one orbit. The reflections group W contains three rotations 
of angles 2rmr/3,m £ {1,2,3} and three reflections C B z ^g2im7r/3^ 
When choosing S = {^1,52} = {e"*'^/^, e*'^/^}, C is a wedge of angle tt/3 in 
the positive quadrant of the plane and 

U«,6Vk{0 < {wsi,x) <e,x£ wC} 

is the hachured part of the following picture 




With some manipulations, one writes (1/0(6), si) as 



L{(eW3,0,)>O} 



dLO((e, si)) + r l|^,..,e.>>o}di°((e, S2)) 
Jo 



+ 



4. A MULTIVOQUE SDE 

Maximal monotone multivoque operators were shown to be a quite effi- 
cient tool to deal with SDEs having singular drifts (see references in jlj). 
They also opened the way to study the strong existence and uniqueness of 
diffusions that are normally reflected on the boundary of their state spaces 
(multivoque Skorohod problem). For instance, one deduces from results 
derived in [Ij that there exists a unique strong solution of the following 
n-dimensional SDE: 



(5) 



Xt = Bt + n{Xt)dkt{X) 



8 



REFLECTED BROWNIAN MOTION IN WEYL CHAMBERS 



where X is a continuous process valued in some convex closed domain D C 
M", B \s a. n-dimensional Brownian motion, n{x) belongs to the inward 
vector field at x (see below) and k \s a, continuous process of finite variation 
such that 

/ ^{Xs(iD}d\k\s = ^- 
Jo 

In this spirit, a SDE similar to ([5]) may be derived for 7r{Q) however with 
oblique reflections. To proceed, recall from [T] that the field of unitary 
inward normal vectors n{x) to C at x € dC is defined by 

{{n{x),x — a), for alia G C}. 

Proposition 2. For any hitting time of the boundary t, dL^ is the sum of 

an inward normal vector at iT{Qt) and an inward normal vector at Oy. 

Corollary 2. The RBM in Weyl chambers satisfies 

TT{Qt) = vr(eo) +Bt+ f[n{TT{Qs)) + no(s)]d^. 

where {lt)t>o is an increasing process and no(s) is an inward vector at Oy 
at time s. 

Proof, the expansion of ^"(0) in the dual basis transforms ([2]) to 

n 

7T{et) = TT{eo) + Bt + ^Yi^i. 

i=l 

Now expand 

n 

and note that any simple root is an inward normal vector at any vactor lying 
in the corresponding hyperplane, therefore it is so at the vertex Oy since Oy 
is the intersection of the walls of the simplex. Moreover, 7r(B) cannot visit 
two or more than two hyperplanes at the same time since does the same 
([6], p. 255). Therefore there exists one and only one index 1 < i < n such 
that (7r(Gf),Sj) = at any hitting time of the boundary so that only Yi 
survives in the drift at that time. As a matter of fact, it only remains to 
prove 

Lemma 1. > for any i ^ j. 

Proof of the Lemma: this is a tricky linear algebra exercise and mainly 
follows from the fact that {sj,Si) < for any i ^ j (|7j p. 9). In fact, 
let A = {{si,Sj))ij = AA^ be the Gram matrix associated with S and let 
= be its inverse (A is invertible since S" is a basis). Then, 

since A is symmetric (hence may be diagonalized in an orthogonal basis) and 
its eigenvalues are strictly positive (it is a positive invertible matrix since 
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A is invertible), the following functional calculus identity makes perfectly 
sense: 



Since —A has positive offdiagonal entries and negative diagonal ones, then 
e"*"^ has positive entries since e~^^ = eK-^+<^^)Q-<^i where c = maxi<j<„(si, Sj 
The Lemma is then proved. ■ 
The Proposition easily follows and the statement of its Corollary is obvious 
with k = Yi{t),n{Tr{Qt)) = {Ci,^i)si and 

at time t such that (7r(0t), Si) = 0. ■ 

Remark. The presence of an inward normal vector at hints to the tran- 
sience to infinity of tt{Q). Indeed Q and 7r(0) has equal length which is a 
Bessel process of dimension > 2, the latter process being transient (|10jj. 
In particular, when R = B2 P-4^)> then C is a wedge of angle 7r/4 and 
according to the convention adapted in [13j on the reflections angles (see 
p. 409), these are negative: 9i = 62 = — 7r/4 so that the RBM is transient to 
infinity (see [15] p. 762). 



5. Further developments 
5.1. Particles system. By the occupation density's formula ([lOj), one has 

Yi = linr V / l{o<(s.,«,*e,><e,e.et«c}'^s, 1 < i < n. 

Thus, since tt{x) = w*{x) when x E wC, then 

Yi = — z — hm - / l{o<(s„7r(0,)><e}as 



2 e^o e 







and again the occupation density formula yields 

Yi = ^ '1 lim - / l{o<(si,7r(es)><e}^'S- 



2 e^o e 







Since (sj, vr(G)) is a continuous semimartingale with brakets process ((sj, Sj)t)j , 
then (Corollary 1.9, Ch.VI in [TO]) 

L0((s,,7r(e))) = 2y, 
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where L^((sj, 7r(G)) is the local time at of the real-valued semimartingale 
(si,7r(0)). Hence 

n 

7r(eo = vr(eo) + 5^ + ^y^e^ 

1 " 

= 7r(eo)+i?i + -j;LO((s,,7r(G)))6. 

i=l 

The lastly displayed SDE resembles 

n 

dXt = dBt + Y,Ll{{n^,X))V, 

i=l 

considered in [11] while studying a system of Brownian particles that reflect 
on hyperplanes orthogonal to n-dimensional unit vectors (nj)j. Above, Vi 
lies in the hyperplane nf- for all 1 < z < n so that the framework presented in 
[TT| is slightly different from ours since the RBM in Weyl chambers performs 
its (oblique) reflections along the dual basis vectors. 

5.2. Multidimensional skew Brownian motion. The skew Brownian 
motion may be defined as a reflected Brownian-like process whose excursions 
from change sign according to independent Bernoulli variables of the same 
parameter < p < 1 (see ^Hj p. 434). It is the unique strong solution of ([9] 

p.43i)g 

(6) Xt = Xo + Bt+^^^^^L^,{X) 

so that one recovers the RBM for p = 1. Quite complicated ways were pro- 
posed to define a multidimensional version of the skew Brownian motion ([9] 
p. 455). Here, one may construct a continuous process valued in V based on 
the excursions of tt{Q) from the hyperplanes delimiting C. Easily speaking, 
the skew process has continuous paths, evolves like 7r(B) in the interior of 
chambers wC, w £ W and either reflects or transmits at hyperplanes. A 
stochastic analysis of this process will be the topic of a subsequent paper. 

Acknowledgments: the author wants to thank Professor D. Bakry for 
his assistance in proving Lemma [H He is also grateful to his colleague J. 
Zender for providing him with the Latex code used to draw the picture. 
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